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Abstract 


This  report  discusses  the  construction  of  a  tandem  mirror  plasma  propulsion  facility, 
the  numerical  modelling  of  the  hybrid  plume  exhaust,  and  rf  heating  of  the  plasma.  A 
preliminary  experiment  of  the  ICRH  (Ion  Cyclotron  Resonance  Heating)  heating  of  plasma 
ions  was  carried  out.  For  2.0  kW  ECRH  (Electron  Cyclotron  Resonance  Heating)  power 
injected  into  the  central  cell  and  10  kW  ICRH  power  into  the  end  cell,  the  results  obtained 
from  the  probe  in  the  central  cell  are:  ne  —  2.5  x  1016  m~3  and  Te  =  80  eV  (928,000  K) 
in  the  central  cell.  The  estimated  values  in  the  end  cell  are:  nj  =  1.25  x  1017  m~3  and 
Tj  =  500  eV  (5,797,000  K).  The  power  conversion  efficiency  was  about  80%.  The  results 
from  time  dependent  3-D  three  fluid  numerical  modeling  indicate  that  a  boundary  layer 
can  be  formed.  The  formation  of  this  layer  is  strongly  dependent  on  neutral  jet  geometry 
and  injection  angle.  The  ICRH  heating  of  plasma  was  modeled  numerically  and  power 
absorption  efficiency  is  about  50%.  Analytical  analyses  was  done  on  slab  geometry. 

1  Concept 

A  hybrid  plume  is  one  in  which  the  exhaust  fluid  is  a  stratified  mixture  of  hot  plasma 
and  neutral  gas  in  a  magnetic  field.  The  macroscopic  properties  near  the  wall  boundary 
(i.e.,  temperature,  density,  fluid  velocity,  thrust,  etc.)  exhibit  a  drastic  radial  variation 
over  small  distances.  The  exhaust  is  extremely  hot  in  the  core  and  relatively  cold  near  the 
edge  [1,2]. 

Such  a  hybrid  p'ume  can  be  produced  by  surround'1  '  e  hot  plasma  exhausted  from 
the  end  of  a  tandem  mirror  magnetic  confinement  device  with  an  annular  hypersonic  gas 
jet  coaxial  with  the  plasma  or  with  various  pitch  angles  [3,4,5].  The  resulting  plasma-gas 
structure  has  useful  applications  in  rocket  propulsion  of  wide  power  levels  where  I>p  and 
thrust  can  be  continuously  varied. 
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2  Accomplishments 


During  the  past  year  we  have  devoted  our  major  effort  to  the  construction  of  the 
tandem  mirror  magnetic  plasma  confineu.cnt  device.  The  device  was  assembled  in  a  short 
eight-month  period  and  the  facility  became  operational  in  March  1989.  An  overall  view 
of  the  device  is  shown  in  Fig.  1.  A  microwave  created  hydrogen  plasma  has  been  heated 
bv  rf  power  to  a  temperature  for  electrons  of  80  eV  which  translates  to  928,000  K  and  an 
estimated  temperature  for  ions  of  500  eV  or  the  equivalent  of  5,797,000  K. 

After  years  of  wrestling  with  numerical  instabilities,  particle  conservation  and  outflow 
boundary  problems,  and  result  interpretation,  meaningful  and  positive  hybrid  plume  char¬ 
acteristics  have  been  produced  by  computational  modeling.  Preliminary  results  indicate 
that  a  boundary  layer  can  be  found  with  proper  injection  of  neutral  gas.  This  reconfirms 
previous  steady  state  studies  but  with  much  more  interesting  aspects. 

2.1  Experimental  Facility 

The  tandem  mirror  device  was  built  under  the  instrumentation  program  with  addi¬ 
tional  support  from  JPL/NASA.  The  detailed  construction  work  on  the  device  has  been 
described  in  the  final  technical  report.  ft.  is  briefly  summarized  below. 

The  tandem  mirror  plasma  rocket  experimental  device,  a  scaled  down  compact  version 
of  the  fusion  device,  is  very  complicated  in  structure  and  operation.  It  involves  sophisti¬ 
cated  mechanical  components,  ultra-high  and  delicate  vacuum  system,  magnetic  coil  and 
cryogenic  system,  high  current  power  supplies,  microwave  power  transmitters,  high  power 
and  high  voltage  radio  frequency  wave  power  transmitters.  The  involvement  of  diagnos¬ 
tic  equipment  is  even  more  extensive.  A  Langmuir  probe  can  only  take  measurements  in 
the  edge  of  the  plasma  where  the  temperature  is  low;  it  cannot  measure  ion  temperature. 
More  equipment  is  needed,  such  as  a  microwave  interferometer  for  measuring  electron  den¬ 
sity,  from  which  ion  density  can  be  inferred  assuming  charge  neutrality;  a  laser  Thomson 
scattering  system  or  a  visible  light  spectrometer  for  measuring  ion  temperature;  a  particle 
beam  lor  measuring  plasma  potential;  a  laser  fluorescence  system  for  measuring  neutral 
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density,  a  charge  exchange  analyzer,  retarding  field  energy  analyzers,  radiation  detectors, 
diagmagnetic  ioops,  etc.  In  terms  of  operation,  the  pulsing  of  the  magnet,  rf,  mieiowave, 
gas  flow  and  the  diagnostic  systems  have  to  he  synchronized  in  a  proper  time  sequence. 
Therefore  one  of  the  purposes  of  this  research  is  to  develop  methods  of  simplifying  this 
process.  However,  every  diagnostic  is  needed  to  understand  the  plasma  property  fully 
before  this  can  happen.  Because  of  budget  constraints  and  manpower  limitations,  it  will 
take  a  long  time  to  gradually  build  up  this  capability. 

The  magnetic  coil  system  and  vacuum  system  were  assembled  in  August  1988.  High 
vacuum  was  achieved  immediately.  The  magnetic  system  was  operated  at  low  power  and 
microwave  discharge  was  initiated.  To  reach  this  point  a  lot  of  hard  work  was  involved,  such 
as  extensive  metal  polishing  of  the  coil  casing  and  vacuum  chamber,  repeated  electrical 
testing  of  the  coil,  vacuum  testing  of  the  chamher  and  coil  casing,  and  careful  assembly, 
alignment  and  calibration  of  the  system. 

From  September  1988  through  January  1989  we  worked  to  complete  the  construction 
of  the  facility.  To  achieve  that  the  following  tasks  were  done:  In  our  laboratory  we 
constructed  a  power  control  unit  which  was  installed  and  tested.  We  installed  the  high 
current  cables.  The  cryogenic  system  was  installed  and  tested  and  a  resistance  drop  of 
factor  5  was  achieved  when  the  coil  was  cooled  down  to  liquid  nitrogen  temperature.  Two 
rf  transmitter  units  were  brought  up  to  operation.  A  laser  diagnostic  system  was  set  up 
and  optical  alignment  was  tested.  The  end  chamber  was  modified  by  installing  two  large 
rectangular  ports  and  small  ports  for  the  implementation  of  the  laser  fluorescence  system. 
A  second  pumping  line  was  installed. 

Testing  of  the  system  began  in  February  1989.  A  giant  step  was  taken  in  energizing 
the  magnetic  coils  with  high  current.  Because  of  the  urgent  need  of  bringing  the  system 
to  full  operation  rapidly,  we  took  the  risk  of  bringing  the  current  to  about  60%  of  the 
designed  value  quickly  causing  one  coil  to  short.  It  had  to  be  repaired.  The  mirror  coil 
set  and  the  set  of  central  cell  coils  and  booster  coils  are  powered  by  two  separate  power 
supplies.  The  power  supplies  were  found  to  interfere  with  each  other  as  shown  by  the 
current  waveform  in  h  ig.  2a.  1  he  power  supplies  had  to  be  adjusted  to  obtain  the  current 
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waveform  in  Fig.  2b. 

The  next  step  was  to  inject  rf  power  into  the  chamber.  It  was  found  that  the  microwave 
and  rf  powers  could  not  be  operated  simultaneously.  The  microwave  is  in  resonance  with 
electron  cyclotion  frequency  at  low  magnetic  field  and  rf  is  in  resonance  with  ion  cycle  tron 
frequency  at  high  field.  The  microwave  antenna  had  to  be  relocated  to  the  central  cell 
where  the  field  is  a  factor  of  10  lower  than  the  field  in  the  end  cell.  The  end  cell  antenna 
was  installed  in  the  center  plane  of  the  cell  which  blocked  all  ports  and  had  to  be  modified 
and  reinstalled.  Accomplishing  these  two  modifications  was  a  major  operation  because 
the  coils  had  to  be  separated.  Again  due  to  the  race  with  time,  the  rf  power  was  injected 
before  the  circuit  was  optimally  tuned.  The  coupler  of  the  transmission  line  was  damaged 
by  arcing  from  the  reflected  power.  After  the  system  was  checked  out  and  repaired,  we 
spent  weeks  tuning  the  rf  matching  circuit.  In  the  meantime  a  semi-automatic  control 
system  was  installed.  The  pulse  length  of  charging  the  magnet  can  be  limited  bv  a  time 
sequencer  to  eliminate  accidental  human  error. 

The  pressure  in  the  chamber  was  found  to  be  rising  too  fast  during  the  test.  This 
indicated  that  an  ultra-high  vacuum  procedure  of  baking  and  glow-discharge  cleaning  had 
to  be  taken.  This  was  a  long  and  tedious  process  which  was  avoided  again  in  favor  of 
speeding  up  the  program.  Also,  all  the  mechanical  pumps  had  to  be  serviced.  All  the 
diagnostics  were  taken  down  in  order  to  carry  out  thorough  checking  and  testing  with 
repetative  baking  and  cleaning  to  establish  a  base  vacuum  condition. 

The  experience  of  building  such  a  compact  fusion  type  device  might  be  proven  useful 
in  the  future  for  space  exploration.  It  is  inheritantly  complex  and  costly.  To  design, 
construct  and  make  it  operational  was  quite  a  challenging  task  and  it  is  in  fact  a  major 
part  of  scientific  undertaking. 


2.2  Initial  Experiment 


I  he  plasma  heating  experiment  has  been  carried  out  during  the  testing  period  by 
back  tilling  the  chamber  to  a  pressure  from  5  x  HR3  Torr  to  20  x  10~3  Torr  at  various 
field  intensities,  rf  power  levels  and  pulse  lengths.  The  only  diagnostic  we  had  was  the 
Langmuir  probe.  The  Langmuir  probe  is  a  fundamental  diagnostics  and  provides  a  basic 
and  reliable  measurement.  However,  it  has  a  severe  limitation  and  fundamental  difficulties. 
One  has  to  be  very  careful  to  prove  that  the  measurement  is  correct.  The  fundamental 
limitation  is  that  the  probe  comes  in  contact  with  the  plasma  and  it  can  be  burned  easily. 
The  probe  can  also  pick  up  any  electrical  noise,  fluctuations  and  is  particularly  susceptible 
to  the  interference  of  rf  power.  The  result  could  be  very  deceiving.  For  a  long  time  we 
were  skeptical  on  the  probe  measurements.  Such  prudence  was  proved  to  be  correct  as 
discussed  later. 

Figure  3a  shows  the  probe  measurement  of  an  ECRH  discharge.  The  top  trace  is  the 
waveform  of  microwave  power  injected  into  the  chamber.  The  middle  trace  is  the  probe 
sweeping  voltage  waveform.  The  bottom  trace  is  the  probe  reading  of  the  plasma.  During 
each  sweep  the  probe  signal  was  found  to  be  modulated  by  the  ECRH  power  which  turned 
out  to  be  a  pulsed  system.  The  plasma  was  not  sustained.  A  constant  source  of  microwave 
power  is  tnereore  needed.  We  wfill  modify  the  unit  in  the  future  because  we  cannot  afford 
to  take  the  risk  of  shutting  down  the  experiment  at  the  present  time. 

As  shown  in  Fig.  3b  there  is  a  very  strong  probe  signal  when  rf  power  was  injected. 
The  signal  is  about  a  factor  of  10  larger  than  the  microwave  discharge  which  is  shown  by 
the  small  ripples  on  the  same  baseline.  Although  we  later  found  that  this  is  much  larger 
than  the  pickup  without  plasma,  we  had  doubts  that  the  signal  was  real.  A  long  study 
was  carried  out  on  the  structure  of  the  probe,  electrical  isolation  and  shielding,  ground 
loop  and  filtering  until  a  signal  showing  the  current-voltage  characteristics  typical  of  a 
Langmuir  probe  was  obtained.  Such  a  signal  is  shown  by  the  second  trace  in  Fig.  4.  This 
gave  us  confidence  in  the  measurement.  I  he  electron  temperature  can  be  deduced  from 


F  ;g  .5.  Langmuir  probe  measurements  of  plasma:  (a)  with  microwave  discharge;  (b) 
with  low  power  rf  heating  at  10  ms  pulse  width.  Top  and  middle  traces  show 
microwave  and  probe  sweeping  voltages. 
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The  Charaterist ics  of  Langmuir  probe  shown  by  the 
middle  trace.  The  top  trace  is  the  microve  pulse 
and  the  bot  o'u  trace  is  probe  sweeping  voltage. 
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the  following  equation  6j 

/  dl 

Te  =  c(I  —  I,x)  j  —  (1) 

where  I9X  is  the  ion  saturation  current.  The  electron  density  is  given  approximately  by 


I  ax  j  Wle 

n--Tc\i\fT) 


(2) 


where  .4  is  the  cross  sectional  area  of  the  probe  tip  which  is  a  sphere  of  2  mm  in  diameter, 
7',  is  the  ion  temperature,  and  e  is  the  charge  and  mass  of  electrons. 


2.3  Exhaust  Modeling 

The  exhaust  modeling  has  been  carried  out  by  solving  steady  state  MHD  equations.  A 
time  dependent,  three  fluid  (electron,  ion  and  neutral)  code  by  solving  transport  equations 
in  cylindrical  geometry  including  interactions  between  species  has  been  developed.  We 
began  to  learn  the  dynamic  behavior  of  the  plasma-gas  system.  The  theory  and  the 
simulation  method  are  described  below. 

The  macroscopic  quantities  (n,V,  and  T)  are  modeled  via  a  fluid  representation.  The 
basic  equations  can  be  obtained  from  the  work  of  Braginskii  [7].  Braginskii  assumed 
that  there  was  no  processes  that  convert  particles  of  one  species  to  another,  so  the  right 
hand  side  of  his  particle  conservation  equation,  for  example,  is  set  to  zero.  We  have  such 
processes,  and  must  generalize  the  equations  appropriately.  We  need  fluid  equations  for  the 
conservation  of  particles,  momentum,  and  energy.  The  equations  to  be  presented  will  be 
for  the  plasma  ion  species.  The  equations  for  electron  species  are  identical.  The  equations 
for  the  neutral  gas  are  the  same,  except  for  the  momentum  conservation  equation,  which 
lacks  the  Lorentz  force  term.  To  close  the  fluid  equations,  we  need  an  equation  of  state  to 
specify  the  relation  between  pressure  and  temperature. 
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The  Transport  Eanoti^ns 

The  equation  for  conservation  of  particles  is  (Eq.  (1.11)  of  Braginskii  in  reference  7) 

-  V  ■  (nV)  =  Sp  .  (3) 

Here  n  is  the  particle  density  and  V  is  the  mean  fluid  velocity,  defined  as 

V  =  (v)  .  (4) 

v  is  the  exact  particle  velocity,  and  the  ( • )  denotes  an  average  over  velocity  space.  Sp  is 
the  particle  source  term,  to  be  described  in  detail  later. 

The  equation  for  the  conservation  of  momentum  is  (Eq.  (1.12)  of  Braginskii  in  reference 

’i’) 

~  (mnV)  4-  V  •  (mn  (vv))  —  en  -i — V  x  =  R  .  (5) 

Here  v  =  v1  V,  the  exact  particle  velocity,  and  v'  is  a  random  velocity. 

Ignoring  viscosity,  this  becomes: 

EtVxBj  =  R  .  (6) 

where  m  is  the  mass  of  a  particle,  c  is  the  speed  of  light,  p  is  pressure,  E  is  the  electric 
field  and  B  is  the  magnetic  field.  R  is  the  source  of  momentum,  to  be  described  in  detail 
later.  Note  that  both  of  these  equations  are  in  conservative  form,  and  this  is  the  way  that 
they  are  numerically  differenced  in  the  code. 

For  conservation  of  energy,  we  start  from  the  heat-balance  equation,  (Eq.  (1.23)  of 
Braginskii  in  reference  7) 

3  PnT  /  3  \ 

2  -  V-  [~nTVj  +nTV-V  +  U:  W  +  V.q  =  Q  (7) 

Here  q  is  the  heat  flux  and  n  is  the  viscosity  tensor.  Q  is  the  heat  source  term,  to  be 
described  in  detail  later.  Ignoring  viscosity,  this  becomes 

~  V'  {1UTV)  ^nTV'  V  +  =  Q  <8> 


( mnV )  -  V  ■  (mnV  V)  4-  Vp  -  en 
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The  second  term  can  be  divided  as  follows: 

TV  -nV  +  nV  -VT  . 

The  first  term  above  can  be  combined  with  the  appropriate  time  derivative  term  of  Eq.  (8) 
to  give  the  continuity  equation  scaled  by  T. 

?r(|,V„v)=?TSp  ,9) 


Writing 


Q-=Q~-TSP. 


and  dividing  Eq.  (8)  by  (3/2)  n  leaves  us 


^  +  V-VT+ ?TY-V  =  .  (11) 

The  second  and  third  terms  can  be  rewritten  as 

V  -  (TV)  -  -TV  ■  V  .  (12) 

3 

This  leads  to  the  following  equation  for  temperature: 

^  -  v.(rv)- bv  v=  v.q)  .  (13) 

Heat  conduction  is  dominated  by  the  electrons,  and  thus  q  is  only  present  in  the 
electron  temperature  equation,  or,  q  =  qe.  The  form  for  qe  is: 

qe  =  Ke±Vx(kTe)  .  (14) 

Here  k  is  Boltzman’s  constant.  Ke  is  the  perpendicular  electron  thermal  conductivity  and 


is  given  as 


Ke  =4.7 


meu*eTe 


Here  ujce  is  the  electron  gyrofrequency  and  re  is  the  basic  collisional  time. 
To  close  the  fluid  equations,  the  equation  of  state  for  pressure  is 


p  =  nkT  . 
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The  Source  Terms 


The  source  terms  are  from  various  collision  processes  with  the  neutral  gas,  including 
atomic  reactions.  The  details  of  the  various  source  terms  will  now  be  presented. 

The  atomic  reactions  included  here  are  ionization,  charge-exchange  and  recombination 
and  are  characterized  as  follows. 

Electron  Impact  Ionization: 

e  +  H°  —*  H+  -+-  2e  ,  (17) 

Resonant  charge-exchange: 

Hot  cold  cold  hot 

H+  +  H°  — +  H+  +  H°  ,  (18) 

Collisional-radiative  recombination: 

e  +  H+  —  Hm  +  hv  ,  (19) 

and  electron  three-body  recombination: 

e  +  H++e-*H+e.  (20) 

In  the  recombination  reactions  above,  the  names  describe  the  way  the  energy  released 
by  the  recombination  is  carried  away.  H *  is  an  excited  state,  and  hv  represents  photons. 

The  particle  source  term,  Sp ,  consists  of  sources  from  atomic  reactions  between  the 
plasma  and  the  neutral  gas,  which  can  be  broken  into  three  terms 

Sp  —  S\on  +  Sex  +  Srec  .  (21) 

Here  5ion  is  the  source  of  particles  from  electron  impact  ionization,  5CX  is  the  source 
of  particles  from  resonant  charge-exchange  and  Srec  is  the  source  of  particles  from  the 
recombination  reactions. 


The  specific  particle  source  terms  have  the  following  form: 

nin2(crv)rx  .  (22) 

nl  and  n2  are  determined  by  the  species  involved  in  the  reaction.  (crv)rx  is  the  reaction-rate 
for  the  specific  reaction.  The  reaction-rates  have  been  previously  calculated  and  stored  in  a 
lookup  table  to  be  used  by  the  code.  The  cross-section,  cx,  is  temperature  dependent.  The 
average  reaction-rates  (crv)  was  calculated  assuming  a  Maxwellian  velocity  distribution  for 
both  species  involved  in  the  reaction. 

The  momentum  source  term  has  the  same  form  as  the  particle  source  term.  Specifi¬ 
cally, 

R  =  /?ion  +  RcX  +  Rrec  •  (23) 

The  R's  represent  the  source  of  momentum  from  the  various  atomic  reactions, or: 

Rrx  =  (mVjj  nxn7{crv)rx  .  (24) 

In  addition,  the  R  term  contains  elastic  collisions  between  the  plasma  and  neutrals, 
i.e.,  collisions  that  do  not  lead  to  atomic  reactions.  These  elastic  collisions  are  commonly 
referred  to  as  friction.  The  frictional  momentum  transfer  can  be  modeled  as  (for  ions 
scattering  off  of  neutrals) 

R{in  =mn(K,-  Vi)  u>in  ,  (25) 

where  U{n  is  the  collision  frequency  for  scattering  of  ions  by  neutrals. 

From  momentum  conservation  we  know  that  Rft\n  =  Electrons  also  scatter 

elastically  from  the  neutrals,  but  in  the  hybrid  model,  this  is  taken  into  account  in  the 
term  for  the  plasma’s  resistivity  (to  be  described  later). 

The  last  source  term  to  consider  is  the  source  term  for  internal  energy,  Q.  This 
includes  the  internal  energy  brought  in  via  atomic  reactions, 

Q  =  Qion  +  Qcx  +  Qrec  ■  (26) 
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Each  Q  has  the  form 


Qrx  =  (^nkT^j  nin2(<Tu)rx  . 


(27) 


For  the  ions  there  is  an  additional  term  for  ion  heating  from  the  electrons 

3 me  nk 


Qt  = 


m,  re 


—  (Te  -  T, )  . 


(28) 


Electromagnetic  Terms 

All  that  is  left  is  to  find  equations  for  E  and  B.  In  a  hybrid  plasma  model,  both  come 
from  the  electron  momentum  equation. 

In  the  limit  of  small  electron  inertia,  the  electron  momentum  equations  reduced  to 

E  = - - - ue  x  B  +  tjJ  ,  (29) 

en  c 

where  rj  is  the  plasma  resistivity.  The  last  term  of  this  equation,  77J,  is  simply  the  R 
term  in  the  notation  used  for  the  ion  momentum  equation.  In  other  words,  the  friction  of 
electrons  scattering  off  of  ions.  As  mentioned  above,  this  term  must  also  account  for  the 
scattering  of  electron  off  of  neutrals.  This  is  taken  into  account  in  the  r/  term. 


*ltotal  —  *7en  "t"  Vei  •  (30) 

The  electron  momentum  equation  can  be  considered  the  Ohm’s  Law  for  the  system.  It 
should  be  emphasized  that  the  velocity  in  the  Lorentz  term  of  this  equation  is  the  electron 
velocity,  ue,  not  the  fluid  velocity  V.  These  two  velocities  are  related  through  the  current, 

J  =  en  (V  -  ue )  ,  (31) 

or, 

ue  =  -—  J  +  V.  (32'' 

en 

Thus,  the  ue  x  B  term  can  be  expanded  to  the  form 

V  x  B  -  —  J  x  B  (33) 

en 
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This  second  term  is  commonly  referred  to  as  the  Hall  current  term.  In  a  more  familiar 
form,  the  electron  momentum  equation  becomes 


E  +  VxB  =  i)Jh - (JxB  —  Vpe)  . 

en 


(34) 


However,  we  will  use  Eq.  (29)  and  Maxwell’s  equations  to  transform  the  electron  momen¬ 
tum  equation  into  an  equation  for  B. 

The  radiation- free  limit  of  Ampere’s  law  is 


V  x  B  =  —  J 
c 


(35) 


Faraday’s  law  is 

=  ,36) 

Introducing  the  magnetic  vector  potential  (Vx  A  =  B)  in  the  Coulomb  gauge,  and  utilizing 
Maxwell’s  equations  we  can  transform  the  electron  momentum  equation  to 

A  =  +  ue  x  B  -  —  q(V  x  B)  .  (37) 

en  47r 


The  9  component  of  this  equation  will  give  us  both  Br  and  Bz.  Finally,  if  we  take  the  curl 
of  the  electron  momentum  equation,  and  again  utilizing  Maxwell’s  equations,  we  obtain 
the  equation  which  was  used  to  advance  Bg  numerically. 

The  electron  velocity  components  are  given  via  the  relationship  between  the  fluid 
velocity  and  current.  Specifically, 


Utr 

Uer 


—  k,r  + 

=  14.+ 

=  viM  - 


dBg 


Airen  dz  ' 
-V2A' 


e , 


4nen 
c  1  d(rBg) 
47ren  r  dr 


(38) 


This  completes  the  description  of  the  equations  for  numerical  modeling  of  the  system. 
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Boundary  Conditions 


To  complete  the  formulation  of  the  problem  one  must  specify  consistent  boundary 
conditions  (EC's)  for  the  equations  discussed  above.  The  geometry  is  a  cylindrically 
symmetric  duct.  EC's  must  be  specified  at  the  inflow,  zmtn,  outflow,  zmax,  r  =  0,  and  the 
wall,  rmax. 

The  inflow  is  assumed  completely  specified  for  the  plasma  and  held  fixed  in  time.  The 
same  is  assumed  for  the  neutral  gas  for  axial  injection.  If  radial  injection  is  used,  zmin 
becomes  an  outflow  boundary  for  neutrals,  also  to  be  described  later. 

At  the  origin,  r  =  0,  all  quantities  axe  assumed  symmetric  in  0  which  leads  to  the 
BC  d/dr  —  0.  This  holds  for  all  quantities  except  Aq.  .The  BC  for  A$  is  Dirichlet 
zero.  L’Hopital’s  rule  is  utilized  to  ensure  the  fields  derived  from  the  vector  potential  is 
appropriate. 

At  the  wall,  rmax,  Aq  is  specified.  For  the  fluid  quantities,  most  axe  set  as  Neumann 
zero.  The  exception  is  temperature.  In  some  cases  the  temperature  is  specified. 

The  most  difficult  BC' s  occur  at  zmax.  Conventional  outflow  boundaxy  conditions 
are  used  for  the  neutral  gas,  e.  g.,  simple  interpolation  back  in  z.  When  the  same  BC’s 
are  attempted  for  the  plasma,  the  computation  failed  eventually  due  to  the  problem  that 
the  plasma  fluid  reflects  (numerically)  from  the  outflow  boundary.  Many  different  BC's 
were  attempted  for  the  plasma  fluid  at  zmax.  The  best  results  were  obtained  by  simply 
interpolating  back  up  a  flux  surface  (magnetic  field  line),  and  setting  the  fluid  quantities 
constant  along  the  field  line  just  at  zmax ■  In  this  sense  the  magnetic  field  line  is  approxi¬ 
mating  a  characteristic,  even  though  in  actuality  it  is  not.  To  date,  the  best  results  have 
been  obtained  using  the  above  BC's  for  the  system. 
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Discussion 


Our  problem  is  many  degrees  more  general  them  a  simple  one  fluid  MHD  problem  and 
is  therefore  many  degrees  more  complex.  Instead  of  solving  one  set  of  fluid  equations,  we 
are  solving  three  sets  of  fluid  equations  -  one  set  for  ions  and  one  set  for  electrons  plus 
the  neutral  transport  equation.  Instead  of  steady  state  we  are  solving  the  time  dependent 
problem.  The  reactions  between  species  were  taken  into  account.  The  coordinate  system 
is  cylindrical  instead  of  a  two  dimensional  Cartesian.  The  electron  stream  and  the  ion 
stream  are  governed  by  different  time  scales.  The  plasma  time  scale  is  many  orders  of 
magnitude  faster  than  the  neutral  jet. 

We  were  facing  the  difficulty  of  energy  and  particle  conservation,  numerical  instability 
and  boundary  problems,  particularly  the  outflow  boundary.  If  not  properly  treated  the 
plasma  fluid  will  reflect  numerically  from  the  outflow  boundary  as  discussed  earlier.  A  nor¬ 
mal  characteristic  solution  method  using  CDF  failed  because  of  a  diverging  magnetic  field. 
The  problem  was  finally  resolved  by  using  a  combination  of  ghost  cells  and  a  characteristic 
like  solution  based  on  magnetic  field  lines. 

The  size  of  the  code  is  very  large  and  computation  time  is  very  long.  At  present  we  are 
carrying  out  a  quick  scan  over  a  wide  range  of  parameter  space  to  survey  the  effect.  Some 
cases  will  be  selected  for  a  focused  study  later.  The  geometry  of  the  tandem  mirror  device 
is  shown  schematically  in  Fig.  5a.  The  expanded  exhaust  section  is  shown  in  Fig.  5b  and 
is  marked  by  the  coordinate  axis  r  —  z,  where  r  =  0  is  the  centerline  of  the  device.  Z  —  0  is 
the  entrance  of  the  plasma  into  the  exhaust.  Figure  6  is  the  density  contour  of  plasma  and 
neutral  with  a  coaxial  injection  of  gas  jet.  The  boundary  layer  was  formed  which  reconfirms 
our  early  steady  state  solution.  Figure  6a  is  the  flow  of  plasma  without  neutrals  at  4ps. 
The  plasma  is  seen  fanned  out  along  the  field  lines.  Figures  6b  through  6d  show  that 
plasma  with  coaxial  neuteral  jet  at  2 p,  3 p  and  4 p  respectively.  The  plasma  was  contained 
and  boundary  layer  was  formed.  However,  the  neutral  jet  was  eroded  downstream.  Figure 
7  shows  density  contour  of  plasma  and  neutrals  when  the  neutral  was  injected  radially. 
The  plasma  was  contained  better  but  neutrals  were  eroded  at  the  tip.  Figure  8  shows  the 
coaxial  injection  at  increasing  neutral  density.  It  shows  that  the  plasma  was  contained  to  a 
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(a) 


Exhaust 


ri2ure  5.  Schematic  picture  of  the  Tandem  Mirror  Plasma  Rocket:  (a)  Side  view 
of  the  device  with  exhaust  duct  where  the  shaded  area  represents 
plasma;  (b)  The  expanded  exhaust  duct. 
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[fl: 


igur«  Plasma  Neutral  interaction  in  the  exhaust, (a)  Plasma  expands  as  the 
field  lines  without  neutral  injection,  (b),(c),  and  (d)  Boundary 
layer  was  formed  with  axial  neutral  Injection,  Neutral  gas  was 
gradually  eroded  at  later  time. 


Plaaaa  neutral  interction  in  the  exhauat,  (a)  Plaama  expands  aa  tha 
field  lina  vithouc  neutral  injection.  (b),(c),and  (d)  boundary  layer 
was  formed  when  neutral  vaa  injected  radially. 


smaller  radius  at  nn  =  10  x  1021  m  3  in  Fig.  8d  than  at  the  lower  density  of  nn  =  1  x  1021 

m~3  in  Fig.  8b. 

« 

In  the  radial  injection  case  the  flow  characteristics  was  further  improved  by  increasing 
the  width  of  the  radial  gas  jet.  The  flow  time  has  been  extended  to  23 fj.  which  is  nearly 
steady  state.  The  contour  plots  at  2,  3  and  4ps  are  shown  in  Fig.  9. 

Figures  10  (al)  and  10  (a2)  show  the  average  radial  and  axial  ion  temperature  profiles 
of  the  plasma  without  injection.  Figures  10  (bl)  and  10  (b2)  are  the  profiles  for  the  plasma 
with  axial  injection.  Figures  10  (cl)  and  10  (c2)  are  the  profiles  with  radial  injection. 

Figure  10  (al)  shows  that  the  tail  of  the  temperature  profile  extends  to  7.7  cm.  Figure 
10  (bl)  shows  that  tail  was  pushed  back  to  5.5  cm  by  the  boundary  layer  created  with  a 
neutral  jet.  The  ion  temperature  was  nearly  the  same,  therefore  it  was  not  deteriorated 
by  the  gas  except  at  the  edge  where  the  temperature  was  reduced  as  expected.  There  was 
almost  no  change  in  axial  profile.  These  effects  can  be  seen  from  the  superposition  of  the 
two  radial  profiles  and  axial  profiles  shown  in  Fig.  11  (al)  and  11  (a2). 

The  change  in  radial  and  axial  profiles  for  the  plasma  with  radial  injection  of  wider 
width  was  significant.  The  superposition  of  these  profiles  with  the  profiles  for  no  injection 
is  shown  in  Figs.  11  (bl)  and  11  (b2).  Radially  the  temperature  reduced  further  at  the 
boundary  but  increased  in  the  center  core.  Furthermore,  the  core  temperature  is  uniform 
to  r  =  1.7  cm  and  drops  sharply  to  a  low  temperature  at  r  =  3.7  cm.  This  unexpected 
characteristic  will  improve  the  propulsion  efficiency  which  is  unique  to  the  system.  The 
axial  profile  in  Fig.  11  (b2)  shows  that  the  high  temperature  plasma  flows  forward. 

The  comparison  of  radial  and  axial  density  profiles  without  gas  injection  and  with  axial 
and  radial  gas  injections  are  shown  in  Fig.  12.  The  changes  are  much  more  appreciable 
than  temperature  profiles.  The  density  at  the  edge  for  axial  injection  increased  as  shown  by 
Fig.  12  (bl)  where  the  temperature  decreased  as  shown  in  F;'T.  11  (bl).  Again  there  is  no 
noticable  change  in  the  axial  profile.  For  radial  injection  the  density  increased  drastically 
both  at  the  edge  and  in  the  downstream  axial  direction  as  shown  in  Figs.  12  (cl)  and  (c2), 
respectively.  The  density  increases  were  about  eight  times  larger  than  in  the  axial  injection. 
Because  of  a  large  difference  in  scale,  presently  these  figures  cannot  be  superimposed. 
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Figure  10.  Averaged  radial  and  axial  ion  temperature  profiles  :  (al)  and(a2) 
without  injection;  (bl)  and  (b2)  with  axial  injection:  (cl)  and 
(c2)  with  radial  injection 


Figure  12.  Averaged  radial  and  axial  ion  density  profies:  (al)  and  (a2)  without 
injection;  (bl)  and  (b2)  with  aiial  injection;  (cl)  and  (c2)  with 
radial  injection 
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These  analyses  show  that  injection  angle,  gas  flow  rate  and  density  are  very  important 
and  need  to  be  studied  thoroughly.  There  is  a  wide  range  of  parameter  space  that  has  to 
be  covered  in  order  to  understand  the  dynamic  properties  of  this  plasma-gas  system. 

This  code  can  be  modified  for  MPD  and  arcjet  studies.  The  theory  described  is 
fundamental  for  any  plasma  propulsion.  Currently  the  theory  used  in  the  code  developed 
for  MPD  and  arcjet  is  the  reduced  one-fluid  MHD  formulation.  In  reality  ail  ihtae  ay  st  :ns 
should  consist  of  at  least  three  fluids  as  discussed  here. 

In  any  plasma  device  the  electron  and  ion  temperature  are  different.  Particularly,  in 
a  device  like  MPD  or  Arcjet,  where  the  plasma  is  produced  through  electrical  discharge, 
the  gas  may  not  be  fully  ionized  and  there  will  be  large  amounts  of  neutral  gas.  The 
degree  of  ionization  and  the  interaction  between  the  species  (eletron,  ion  and  neutral) 
should  be  important.  The  electrons  have  very  high  mobility  and  determine  the  current 
flow.  In  principle  the  thrust  of  these  devices  is  from  the  j  x  B  force.  Therefore  it  is  to  the 
best  advantage  of  these  system  that  the  plasma  behaves  as  a  single  fluid.  Otherwise  the 
thrust  will  be  determined  bv  the  temperature  of  the  ions  which  axe  the  momentum  carrying 
species.  Therefore  the  detailed  study  of  the  microscopic  properties  of  the  plasma  using  the 
time  dependent  and  multi-species  transport  code  may  allow  us  to  learn  the  details  of  the 
current  flow,  magnetic  field  and  energy  distribution  and  the  density  of  the  species  involved. 
This  may  lead  us  to  learn  conditions  for  single  fluid  operation. 

Because  there  is  no  internal  electrode  in  the  exhaust  of  the  hybrid  plume  system,  there 
is  no  internal  boundary.  The  initial  conditions  of  the  plasma  injected  from  the  tandem 
mirror  device  are  known.  To  adapt  this  code  for  MPD  and  Arcjet  studies,  the  boundary 
conditions  at  the  internal  electrodes  and  the  initial  conditions  have  to  be  worked  out.  Our 
experience  is  that  to  incorporate  the  time  dependent  internal  boundary  conditions  is  not 
trivial.  To  assure  correct  physics  and  numerical  stability  requires  a  considerable  amount 
of  effort. 
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2.4  ICRH  Plasma  Heating 


Present  fusion  devices  use  ICRF  power  and  neutrai  jeams  to  heat  the  plasma.  Since 
the  neutral  beam  is  larger  and  much  more  complex  than  ICRF,  we  chose  to  use  the  ICRF 
heating  method. 

T^e  modeling  of  ICRH  heating  of  plasma  has  been  done  and  reported  in  the  past 
which  provided  the  guideline  for  antenna  design.  Analytical  studies  on  slab  models  was 
also  done.  For  the  completeness  of  discussion  they  are  reviewed  here. 

The  analysis  of  the  problem  requires  numerically  solving  for  the  fields  in  the  vacuum 
and  plasma  regions.  A  code  developed  by  B.  McVey  [8]  was  used  to  solve  this  problem  for 
different  antenna  geometries.  The  theory  used  in  this  code  is  briefly  summarized  below. 

Theory 

The  code  analyzes  Maxwell’s  equations 
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and  solved  in  both  the  vacuum  and  plasma  regions.  The  difference  between  the  vacuum 
and  plasma  region  comes  in  the  current  term  J. 
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* 

J  vacuum  ^ 

(45) 

J plasma  X) 

(46) 

where  e  is  the  dielectric  tensor. 

The  code  solves  these  equations  by  breaking  up  the  cylinder  into  several  concentric 
cylindrical  regions,  solving  for  the  fields  in  each  region,  then  matching  the  boundary  con¬ 
ditions  to  couple  to  the  next  region  with  its  different  parameters  of  density,  temperature. 

Strata  Wave  Matching  Geometry 

The  transverse  fields  in  each  region  can  be  expressed  in  .terms  of  Ez  and  Hz,  which 
are  found  to  be 

VACUUM 

Et  =  EJn(ur)  +  E2Kn(ur)  '  (47) 

Hz  =  HJn(ur)  +  H2Kn(vr )  (48) 

PLASMA 

Hz{r)  =  HzlJn(k±ir)  +  y2Ez2Jn(kj_2r)  (49) 

Ez{r)  =  Ez2Jn(kL2r)  +  zxHzlJn(k^r)  .  (50) 

From  these  fields,  the  time  averaged  power  flow  per  cross-sectional  area  can  be  calculated 

by 

P(w/m2)  =  ^Re[E(f)x  H'(r)}-  ar  .  (51) 

The  inward  power  flow  is  then  found  by 

P(r  =  r0)  =  Y,  I"  Re\E*H*  -  E‘H*\dk‘  ■  (52) 

Z  n  d 

The  power  absorbed  by  the  plasma  is  determined  by  the  plasma  currents  flowing  in  phase 
with  the  electric  field: 

p(w/m3)  =  -  Re{E  -r )  .  (53)  , 
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From  the  circuit  point  of  view,  the  power  transferred  out  of  the  antenna  region  is 


’  dv 


Pc(w) =-\y,J  dk*Jb 


The  irnpedence  is  then  found  by 


Z  =  - 


2Pc(w) 


\H: 


(54) 

(55) 

(56) 


Numerical  Results 

Seven  types  of  antennae  have  been  studied  using  the  ANTENA  computer  code  [8]  and 
axe  shown  in  Fig.  13.  Typical  computed  power  absorption  profiles  by  ions  and  electrons 
are  shown  in  Fig.  14.  The  results  are  summarized  in  Table  1.  The  total  input  power  used 
in  these  calculations  is  100  kW. 

The  current  in  the  coils  necessary  to  provide  100  kW  input  power  varies  from  about 
2.5  to  8.5  kA.  All  coils  are  assumed  to  have  circular  cross  sections  with  a  radius  of  1.0 
cm  located  4  cm  radially  from  the  plasma  boundary.  A  reasonable  limit  for  the  current 
density  is  about  2  kA/cm2,  which  is  exceeded  in  the  full  turn,  the  line  current,  and  the 
half  Nagoya  coils. 


Table  1 

Absorption  for  Various  Antennae 


Coil 

/(kA) 

./(kA/cm2) 

Pt(  kW) 

Fi(%  ) 

Fe(%  ) 

Efficiencyl 

Full  turn 

7.204 

2.293 

45.97 

84.53 

15.48 

45.97 

Saddle  coil 

3.342 

1.064 

49.23 

97.85 

2.15 

49.23 

Nagoya  III 

3.526 

1.122 

49.04 

98.50 

1.50 

49.04 

Aperture 

2.738 

0.872 

28.72 

98.48 

1.52 

28.72 

Line  current 

6.735 

2.144 

54.05 

96.04 

3.96 

54.05 

Half  loop 

5.197 

1.654 

54.90 

95.72 

4.28 

54.90 

Half  Nagoya 

8.276 

2.634 

47.73 

54.14 

45.86 

47.73 
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Typical  absorption  profiles. 
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coil  current 
coil  current  density 
total  power  absorbed  by  ions 
fraction  of  total  power  absorbed  by  ions 
fraction  of  total  power  absorbed  by  electrons 
power  calculated  by  the  antenna  code 

fraction  of  power  input  to  antenna  absorbed  by  the  plasma 

The  results  show  that  the  half  loop  antenna  has  the  highest  efficiency  of  all  the  an¬ 
tennae  considered.  The  ions  in  the  system  absorb  96%  of  this  power  while  the  electrons 
receive  the  balance.  This  is  a  desirable  result  since  energy  transfer  from  the  ions  to  the 
neutrals  is  more  efficient  than  from  the  electrons  to  the  neutrals. 

These  preliminary  results  agree  well  with  experimental  ones  which  show  ICRF  heating 
efficiencies  greater  than  55%  [9].  The  theory  of  ICRF  is  not  well  developed.  It  is  a  very 
complicated  problem  which  is  worked  on  by  all  the  researchers  in  the  fusion  community.  We 
are  only  a  small  effort  toward  the  contribution  to  the  understanding  and  the  improvement 
of  a  heating  mechanism.  VVe  will  continue  to  apply  the  new  development  in  the  fusion 
community  to  this  study. 

Based  on  the  relative  efficiencies  of  the  antennae  shown  by  these  results,  we  have 
chosen  two  half-circle  coils  to  provide  ICRF  heating  at  the  initial  stage.  A  full  turn  loop 
was  not  chosen  because  it  seems  to  require  a  higher  current  and  higher  current  density  to 
provide  the  same  input  power.  Two  half  loops  will  be  used  to  maintain  the  symmetry  of 
the  plasma;  a  single  partial  loop  tends  to  increase  the  plasma  current  on  the  antenna  side. 
The  half  turn  also  has  the  advantage  for  our  application  of  being  light  in  mass  compared 
to  some  of  the  other  coils. 

2.5  Analytical  Modeling  of  ICRH  Plasma  Heating 

The  numerical  model  presented  in  the  previous  section  is  lacking  physical  insight. 
The  analytical  model  based  on  a  slab  geometry  was  done  in  order  to  obtain  the  explicit 
dependence  of  heating  on  the  antenna  geometry.  The  power  coupled  from  an  antenna  into 


/  = 
J  = 
Pt  = 

Fi  = 

Fe  = 

Pa  = 
Efficiency  = 
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a  plasma  may  be  calculated  if  both  the  input  antenna  current  and  the  loading  resistance 
of  the  plasrfia  are  known.  The  power  coupled  into  a  plasma  is  calculated  by 

P=l-Rp\I\2  (57) 

where  Rp  is  the  plasma  loading  resistance,  I  the  applied  antenna  current.  Thus  it  is 
necessary  to  calculate  a  load  resistance  for  the  configuration  desired.  Although  the  physical 
geometry  of  the  plasma  rocket  is  cylindrical,  the  calculations  are  performed  using  a  slab 
geometry  to  produce  a  resulting  Rp  that  can  be  intuitively  understood.  The  geometry  is 
seen  in  Fig.  15.  The  results  presented  are  obtained  from  Shepard  [10]. 

The  analysis  is  as  follows.  The  antenna-plasma  coupling  region  is  separated  into  three 
regions:  Region  I  is  the  vacuum  between  the  chamber  wall  and  the  antenna,  Region  II  is 
the  vacuum  between  the  antenna  and  the  Faraday  Shield,  and  Region  III  is  the  plasma 
region.  The  density  profile  in  Region  III  is  assumed  to  be  a  constant  sloped  line.  The  fields 
are  solved  in  each  of  the  three  regions  and  then  boundarv  conditions  are  used  to  couple 
the  fields  together.  Complete  knowledge  of  the  fields  then  can  be  used  to  find  the  power 
coupled  into  the  plasma  by 

Pa  =  J°°  Ey\antJ-(nz)^  =  l-I2ZA  (58) 

where  Ey  ant  is  the  Electric  field  at  the  antenna  surface,  J*{nz)  the  nz -transformed  antenna 
sheet  current,  and  ZA  the  Specific  Radiation  Impedence  of  the  plasma.  ( Rp  =  Re(Z4)) 
Thus  Rp  can  be  found  from  the  fields  and  a  parametric  dependence  of  the  plasma  loading 
dependence  can  be  obtained. 

The  fields  in  Regions  I  and  II  are  found  by  solving  Maxwell’s  equations  in  free  space. 

V,E=-f  (59) 

v  x  b  -  4  +  ^oJ  (6°) 

c 2  at 

REGION  I. 
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Vacuum 


Region  3rd 
Plasma 


Geometry  of  the  wall,  antenna  and  shield  in  slab 


Defining  c  to  be  1  (where  c  is  the  speed  of  light),  the  vacuum  fields  are  found  to  consist 
of  a  linear  combination  of  Transverse- Electric  (TE)  and  Transverse-Magnetic  (TM)  modes: 


E  =  .4iEta/  +  -^Ete 

(61) 

B  =  TiBrA,/  +  ^Bte 

(62) 

Defining  nx 
TM  MODES: 


the  following  fields  are  obtained  for  Region  I: 


Ex  =  i(n2  +  n2z)  cos(nx() 

Ey  =  r«xny  sin(nx£) 

Ez  -  nxnz  sin (nz£) 

Bx  =  0  (63) 

By  =  inz  cos(n*£) 

Bz  =  —iny  cos(nx£) 


TE  MODES: 

Ex  =  0 

Ey  —  -nz  sin(n*£) 

Ez  —  riy  sin {nx£) 

Bz  =  (nv  +  n2z)sin(nxZ)  (64) 

By  —  inxny  cos (nx£) 

Bz  =  inx  n,  cos(nr£) 


REGION  II. 

The  difference  between  regions  I  and  II  is  the  jump  in  the  B-Field  across  the  antenna 
due  to  the  surface  current  J(nv,n2).  Using  the  boundary  condition  Ez{ 0)  =  0,  one  of  the 
two  constants  .42  can  be  found  in  terms  of  the  other  and  the  resulting  field  at  the  Faraday 
shield  is  : 
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Ex(0)  —  iAiny(riy  +nl)cos(nx6) 

Ey( 0)  =  Axnx(nl  4-  n:)sin(nx6) 

E,(0)  =  0  (65) 

B*(0)  =  -.4in*n,(nJ  +  n])  sin (n,6) 

By( 0)  =  i.4invnz(n^  +  n“)  cos(nl(5) 

B2(0)  =  —iAi(riy  -r  n2xn2z)  cos(nxS)  -  J(nv,nz) 


REGION  III. 

The  method  of  solution  presented  are  from  Lam  et  al.[llj.  The  fields  in  the  plasma 
can  be  found  by  solving  Maxwell’s  equations  with  the  plasma  contribution  appearing  in 
the  current  J .  After  Fourier-Laplace  transforming  and  ignoring  the  initial  condition, 
Maxwell’s  equations  are  reduced  to  : 

VxE  =  tu;B 

VxB  =  -iu/R  ■  E  (66) 


w 


here  K  is  the  cold  plasma  dielectric  tensor: 


Kl  - iKz  0 

K  =  |  iKx  Kl  0 

0  0  Kt  ) 


\ 


where 


h\  =  1-  V  -~ 

‘ — '  i±j  “  — 

a 

K 


U>‘ 


uu 


pa 


ni 


pa 


I 


na^2-ni 

‘-E 


Uu 

JJ 


2 

pa 

2 


(67) 


(68) 


where  a  is  the  species  index. 


39 


Expanding  K\\,  it  is  seen  that  the  dielectric  tensor  K  can  be  separated  into  two  terms, 
one  of  ordetO(l)  and  another  of  order  0{\/t).  (Where  e  —  me/m,i  )  Maxwell’s  equations 
now  become 

V  x  E  =  iu>B 

eVxB=  iS,  •  E  (69) 

where 

S.  =  cS  +  Sqzz  (70) 


Sx  = 


9 

ujuj- 

Pl 


oj2  —  u ;2- 

^  C  l 


Expanding  the  Electric  and  Magnetic  fields  in  orders  of  e, 


E  —  Eq  ”  fEj  B  —  B0  ■+■  tBj  -)-••• 


it  can  be  shown  that  Maxwell  equations  reduce  to  : 


dEy 

dBz 


nySx 


Ey  ~  l 


1  + 


n: 


5,  +  n2 


n  y  S  Y 


5, 


b2 


S'2 

J  x 


>x  +■  ni 


n. 


(71) 


(72) 


(73) 


(74) 
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The  problem  can  be  further  reduced  by  examining  the  case  of  no  propagation  in  the 
y-direction.-In  this  case,  Maxwell’s  equations  are  further  reduced  to 


+  nos,  =  o 

Ci,S 

(75) 

where 

f(i)  -  c  2  ~  s± 

Si  +nj 

(76) 

Approximating  the  function  /(£)  bv  : 

m  =  7(0  -  (nz  - 1) 

(77) 

where  7(0  is  the  density  profile. 

The  solution  of  the  above  differential  equation  is  the  Airy  function: 

Ey(0  =  -43Ai(-i/)  -  A4Bi(-i/) 

(78) 

where 

-  (n2  -  1) 

A/2/3 

v/  =  n°m^ 

eB 2 

^0  =  n(£  =  T) 

(79) 

Using  the  Outgoing  Wave  Boundary  Condition  by  taking  £ 

— ►  00,  it  is  seen  that  an 

outgoing  wave  has  the  wave  property  of  e:A=-"t)  _ 

Using  the  trigonometric 

identity 

elfcl  =  cos(kx)  +  1  sin(fcx) 

(80) 

Asymptotically, 

4.,  ,  ^in(I-3/24-f) 

A,( 

cos(  fi^3'  2  4-  M 

(81) 
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Thus 

Ey  =.43[*Ai(-«/)  +  Bi(-»/)]  (82) 

Using  Taylor  Series  Expansion  !2]  for  small  i/,  the  Electric  and  Magnetic  fields  are 
seen  to  be  : 

Ey( 0)  =  D[(i  +  v'3)  -  (i  -  V3 )F(n]  -  1)] 

Bz(0)  =  D{\  -  iV3)FM  (83) 
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From  this  relationship,  the  parametric  dependence  of  the  loading  resistance  can  be 
seen.  The  results  show  that  for  increasing  FM6,  the  loading  resistance  RA  decreases. 
The  parametric  results  can  be  seen  in  Table  2.  It  is  seen  that  to  obtain  the  maximum 
loading  resistance,  the  distance  from  the  antenna  to  the  wall,  d,  must  be  maximized  but 
the  distance  from  the  wall  to  the  Faraday  shield,  <5,  must  be  minimized.  It  is  also  desirable 
to  have  a  relatively  low  density  plasma  and  a  high  magnetic  field.  The  antenna  width  itself 
must  also  be  kept  as  thin  as  possible  to  maximize  the  plasma  loading  resistance. 

Table  2 

Parametric  Dependences  of  the  Loading  Resistance  RA 


Variable 

Significance  of  Variable 

Result  on  R 

d 

d  — Distance  from  Antenna  to  Wall 

d  T,  Ra  T 

£ 

£  - -  W/2  ~~  width  of  antenna 

a,  Ra  1 

n  o 

no  =  density 

no  T,  Ra  i 

B 

B  =  Magnetic  field 

B  t,  Ra  T 

6 

5=Dist.  from  wall  to  Faraday  Shield 

<5  T,  Ra  1 

Discussion 

A  3-D  cartesian  analysis  of  the  antenna  geometry  was  carried  out  using  the  method 
of  Shepard  [10]  and  Lam  [  1 1  ] .  The  resulting  plasma  loading  resistance  shows  a  parametric 
dependence  on  d,  £,  n0,  B,  and  6  which  are  the  distance  from  the  antenna  to  wall,  width 
of  antenna,  plasma  density,  magnetic  field,  and  distance  from  the  wall  to  the  Faraday 
shield,  respectively.  It  is  seen  that  in  order  to  maximize  the  plasma  loading  resistance,  the 
antenna  must  be  placed  as  close  to  the  Faraday  shield  and  as  far  away  from  the  wall  as 
possible.  Other  parametric  dependences  of  the  loading  resistance  indicates  that  a  thinner 
antenna  along  with  a  lower  density  plasma  and  a  higher  magnetic  field  would  yield  a  higher 
loading  resistance. 
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Future  work  in  this  area  would  extend  the  treatment  from  a  3-D  Cartesian  geometry 
to  a  3-D  cylindrical  geometry.  This  would  better  mock-up  the  IPPROC  and  perhaps  yield 
a  more  accurate  picture  of  the  behavior  of  the  loading  resistance.  Hot  plasma  effects  will 
also  be  mocked  up  so  as  to  better  simulate  the  plasma  inside  the  device. 

3  Propulsion  Projection 


The  exhaust  velocity  of  the  plasma  injected  from  the  mirror  device  is  given  by  thermal 
velocity  [12] 


v  = 


(91) 


where  T  is  the  ion  temperature  in  eV,  m  is  the  hydrogen  ion  mass  in  kG,  and  .4  is  the 
mass  number  for  heavy  atoms.  The  velocity  can  be  rewritten  as 


v  =  0.979  x  104 


m/s  . 


(92) 


Let  us  use  100  kW  of  power  as  a  basis  to  find  the  range  of  feasible  performance  parameters. 
Given  the  power  P  in  k\V  we  obtain  the  mass  flow  rate 


2  p  A 

m  =  — —  =  2.090  x  10-8  —  P  mg/sec 


v‘ 


T 


(93) 


and  the  thrust 

f A 

Thrust  =  mv  =  0.205  J  —  P  N  .  (94) 

The  exhaust  cross-sectional  area  of  this  device  is  S  =  0.785  x  10-2  m2.  One  can  obtain 
the  required  plasma  density 


n=s^  =  1-sfflxl0,V£F  m"‘-  (95) 

The  thrust  and  mass  flow  rate  m  are  plotted  as  solid  curves  as  a  function  of  I,v  in  Fig.  16, 
whereas  the  plasma  temperature  and  density  are  plotted  as  dashed  curves.  The  curves  are 
self-explanatory.  The  thrust  decreases  from  100  N  at  /jp  =  1,000  sec  to  0.6  N  at  Itp  = 
30,000  sec.  The  plasma  temperature  T  increases  from  1  eV  to  1  keV.  The  corresponding 
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density  decreases  from  l  x  1021  m~3  to  1  x  10 17  m-3.  These  plasma  parameters  are 
achievable  in  this  device  but  harder  toward  low  temperature  and  high  density.  This  is 
where  the  hybrid  plume  comes  in. 

Figure  17  shows  the  power  and  mass  flow  rate  as  a  function  of  I,p  (solid  curves) 
at  constant  thrust  of  20  N.  Power  increases  from  100  kW  to  2  MW  and  mass  flow  rate 
decreases  from  2  mg/s  to  0.7  mg/'s  for  ranging  from  1,000  sec  to  30,000  sec.  The 
temperature  range  is  the  same  as  in  Fig.  17  for  the  same  Itp  range.  However,  the  density 
is  two  orders  of  magnitude  lower. 

Figure  18  shows  the  mass  flow  rate  and  power  requirement  as  a  function  of  thrust 
with  constant  I,p  of  3,000  sec.  Therefore  the  temperature  is  fixed  at  10  eV.  The  density 
is  increased  by  two  orders  of  magnitude  than  that  in  Figs.  16  and  17. 

The  present  achieved  conditions  put  us  at  the  high  I,p  (23,000)  sec  and  low  thrust  0.1 
N.  We  are  planning  to  work  downward  to  increase  thrust  and  reduce  These  figures 
indicate  the  range  of  parameters  we  would  like  to  achieve. 

4  Physical  Defflciency 

The  properties  of  the  hybrid  plume  are  governed  by  the  physics  of  high  temperature 
hydrogenic  plasma  interacting  with  hydrogen  gas  jets  in  a  nonuniform  magnetic  field.  The 
physics  of  such  a  system  is  not  well  known. 

Theoretically  the  concept  has  been  verified  by  a  steady  state  fluid  model  in  a  uniform 
magnetic  field  with  simplified  assumptions.  This  was  a  step  to  demonstrate  the  principle  of 
the  hybrid  plume.  Dynamic  flow  characteristics  have  to  be  studied  with  a  time  dependent 
model  in  a  nonuniform  field.  In  a  time  dependent  code  proper  treatment  of  the  outflow 
boundary  becomes  very  important.  Now  we  have  achieved  this  second  step  and  we  begin  to 
understand  the  dynamic  behavior  of  the  plasma- neutral  interaction.  Since  the  parameter 
space  is  very  large,  a  thorough  study  has  to  be  carried  out  as  the  next  step. 

From  the  knowledge  gained  in  this  work  we  found  a  new  defflciency  in  the  theory.  The 
pressure  in  the  transport  equation  is  assumed  to  be  scalar.  In  both  MPD  and  arcjet  and 
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Isp(lOOOsec) 

Fig. 16:  Propulsion  parameters, thrust  and 
mass  flow  as  a  function  of  Isp(solid 
curves,  left  scale),  the  corresponding 
plasma  parameters(dashed  curves, right 
scale)  at  constant  power  P  =  100kW. 
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n(m-3) 


Isp(lOOOsec) 


Fig.  17:  Propulsion  parameters,  thrust  and 
power  (solid  curves,  left  scale)  as  a 
function  of  Isp  and  the  corresponding 
plasma  parameters,  temperature  and 
density  (  dashed  curves,  right  scale) 
at  constant  thrust=20  N 
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Fig. 18:  Propulsion  parameters,  mass  flow 
and  P  (solid  curves,  left  scale)  as  a 
function  of  thrust  and  the  correspon¬ 
ding  plasma  parameters(dashed  curves, 
right  scale)  at  constant  lsp=3000sec 
(T=9.39eV). 
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this  concept,  the  pressure  is  actually  a  tensor  and  should  be  treated  as  such.  Therefore  in 
the  fluid  equation  one  should  use  V  ■  p  instead  of  Vp.  The  viscosity  was  ignored  and  should 
be  incorporated  in  the  future.  The  code  should  be  modified  to  take  this  important  physics 
into  account  and  should  be  continuously  generalized  to  become  as  realistic  as  possible. 

Fxperimentally  heating  the  plasma  to  high  temperature  has  been  proven  by  the  large 
fusion  devices.  However,  the  nature  of  the  operation  is  not  the  same  as  for  propulsion 
purposes.  For  a  fusion  device  one  would  like  to  confine  the  plasma  as  long  as  possible 
and  heat  the  plasma  to  extremely  high  temperatures,  whereas  for  a  propulsion  device, 
one  would  like  to  make  the  plasma  flow  at  modest  temperatures,  but  at  much  higher 
densities.  Some  other  problems  such  as  magnetic  detachment  and  asymmetric  operation 
are  important  for  propulsion,  but  have  never  been  attempted  in  fusion  research.  The  entire 
parameter  space  shown  in  Figs.  16,  17  and  18  need  to  be  mapped  out  by  experimental  data 
with  a  clear  understanding  of  the  plasma  properties.  Although  the  rf  heating  of  plasma  has 
been  studied  by  researchers  in  fusion,  both  theoretically  and  experimentally,  the  physics 
is  still  not  completely  understood. 


5  Objective  and  Benefit 

The  objective  and  benefit  are  closely  related  and  described  in  long  and  near  term 
below. 

The  ultimate  long  term  goal  and  benefit  is  to  develop  a  variable  Isp  and  thrust  rocket 
propulsion  device  with  power  conversion  efficiency  better  than  50%.  The  I,p  can  vary  from 
1,000  sec  up  to  30,000  sec.  The  thrust  is  mainly  restricted  by  the  power. 

* 

The  near  term  goal  is  to  achieve  the  plasma  conditions  (temperature,  density  and  mass 
flow  rate)  for  the  wide  propulsion  performance  characteristics  specified  in  Figs.  16,  17  and 
18  and  to  develop  a  simplified  and  reliable  control  method.  The  near  term  benefits  are 
understanding  the  fundamental  physics  of  heating  the  plasma  with  rf  power.  The  device 
provides  a  test  bed  for  other  advanced  propulsion  concepts  such  as  anti  matter  propulsion, 
for  nozzle  studies,  for  high  temperature  space  material  development  and  for  possible  space 
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plasma  experiments.  Both  the  computational  tool  and  diagnostic  developed  in  this  research 
can  be  used  in  research  in  other  areas  such  as  MPD  and  arcjet.  Therefore  the  long  range 
specific  objectives  can  be  stated  in  terms  of  physics  and  technology  as  follows. 

Physics  Objectives 

1.  To  investigate  the  physics  of  plasma  neutral  interaction  in  thin,  hypersonic  boundary 
layers.  Specifically,  what  is  the  impact  of  cold  (neutral)  gas  injection  and/or  mixing 
into  the  exhausted  plasma  as  it  pertains  to  thrust  augmentation,  thrust  efficiency, 
specific  impulse  variability  and  plasma  magnetic  field  decoupling  as  a  function  of 
various  operating  conditions  and  system  configurations? 

2.  To  investigate  plasma  magnetic  Held  detachment  mechanisms  in  magnetic  nozzles. 
Specifically,  what  are  the  conditions  under  which  effective  decoupling  of  the  exhausted 
plasma  from  the  magnetic  field  lines  will  take  place  such  that  useful  thrust  may  be 
produced0 

3.  To  investigate  the  physics  of  RF  and  microwave-heated  plasmas.  Specifically,  how  to 
best  couple  power  to  the  plasma  anti  under  what  conditions  is  it  more  efficient  to  heat 
ions  over  electrons. 

1.  To  investigate  the  characteristics  .»f  tandem  mirror  configurations  operating  asymmet¬ 
rically.  Specifically,  what  the  specific  impulse  t  l,p  i  of  the  exhausted  plasma  is  for  a 
variety  of  operating  conditions  *>f  plasma  density  and  temperature  and  magnetic  field 

!  or.figurations. 


Long-Term  Technology  Objectives 


1.  To  demonstrate  the  viability  of  the  hybrid  plume  rocket  as  a  plasma  propulsion 
scheme.  Specifically,  to  obtain  a  characterization  of  the  operating  envelope  of  such  a 

de ’.'ice. 

2.  lo  investigate  the  technology  of  magnetic  nozzles.  Specifically,  what  fraction  of  the 
total  kinetic  (thermal)  energy  of  the  generated  plasma  can  be  realized  as  useful  thrust 
for  propulsion? 

3.  lo  investigate  the  operation  of  a  tandem  mirror  device  as  a  variable  I,p  piasma  rocket. 

4.  lo  serve  as  a  pathfinder  toward  the  design  and  construction  of  a  flight  unit  to  be 
tested  in  low  Earth  orbit. 

5.  lo  carry  out  efficiency  studies  under  a  variety  of  system  configurations. 

6  Scope  of  Work  for  February  1990  to  January  1991 

[he  long  term  objective  is  to  achieve  the  plasma  characteristics  for  wide  propulsion 
application  as  shown  in  Figs.  16,  17  and  18. 

As  discussed  in  Section  2.2  the  microwave  feedthrough  is  installed  at  the  right  end 
i  m >rt  h  end  t  of  the  central  cell  and  the  rf  antenna  is  in  the  north  end  cell  as  shown  in  Fig.  19. 
I  he  rf  power  heats  the  plasma  ions  and  creates  an  electrostatic  ambipolar  potential  <i>\  in 
the  north  end.  A  second  antenna  will  be  installed  in  the  central  cell  where  the  bulk  plasma 
can  be  heated  by  a  second  rf  transmitter.  Eventually  a  third  antenna  will  be  constructed 
and  installed  in  the  south  end.  Then  the  third  rf  transmitter  will  be  needed  to  heat  the 
plasma  in  the  south  end  to  create  the  electrostatic  ambipolar  potential  <t> 2-  This  completes 
the  tandem  mirror  configuration.  Many  parts  of  the  third  transmitter  were  removed  for 
repairing  the  other  two.  Extensive  parts  and  work  are  required  to  bring  up  this  transmitter 
to  ■  (peration. 
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At  present  we  have  only  Langmuir  probes  installed  in  the  central  cell  which  can  only 
measure  thq  electron  temperature  at  the  cold  plasma  surface.  We  need  to  measure  the  ion 
temperature  for  two  reasons.  The  first  reason  is  that  the  ions  carry  the  momentum  which 
are  the  main  thrust  producing  elements.  1’he  second  reason  is  that  the  potentials  <f>\  and  d> o 
are  proportional  to  the  ion  temperature.  The  proper  method  for  measuring  temperature  is 
laser  Thomson  scattering  which  is  too  expensive  and  complicated.  To  set  up,  operate  and 
make  measurements  with  Thomson  scattering  requires  at  least  a  full  time  scientist.  We 
are  going  to  try  to  determine  the  ion  temperature  from  Doppler  broadening  of  radiations 
produced  by  helium  using  a  monochromatic  spectrometer.  It  is  not  a  simple  task  even 
to  make  Doppler  broadening  measurements  work.  It  needs  very  high  time  resolution  and 
much  faster  time  scanning  than  the  built-in  scanning  system.  We  are  planning  to  use  a 
rotating  mirror  to  disperse  the  frequency.  At  present  we  are  lucky  to  be  able  to  borrow  a 
spectrometer  for  several  months  which  will  allow  us  to  study  such  a  rotating  mirror  system. 

Two  diamagnetic  loops  were  installed  outside  the  chamber  located  at  the  north  end 
cell  and  the  central  cell.  The  test  results  indicate  that  a  signal  was  shielded  by  the  vacuum 
chamber.  Therefore  the  loops  need  to  be  installed  in  the  interior  of  the  chamber.  To  make 
such  an  installation  is  very  complicated.  The  loop  has  to  be  wound  inside  a  tube  and  be 
vacuum  tight,  and  it  also  has  to  be  protected  from  the  plasma.  The  chamber  has  to  be 
opened  at  many  locations  for  the  installation.  The  diamagnetic  loop  measures  the  ratio  of 
plasma  pressure  to  magnetic  pressure  called  3 

3  =  P/(£2  2p)  .  (96) 

Since  the  magnetic  pressure  B2/2p  is  known,  the  overall  plasma  pressure  can  then  be 
determined 

MS)- 

To  measure  the  energy  of  the  ions  md  electrons  in  the  exhaust  we  need  a  retarding 
field  grid  energy  analyzer.  Such  an  energy  anlyzer  developed  for  fusion  research  is  shown 
in  Fig.  20.  The  electronics  for  driving  the  probe  includes  power  supply,  functional  gener¬ 
ators  and  a  multi-channel  digitizer.  The  probe  itself  is  not  only  expensive  but  also  very 
complicated  and  delicate.  The  probe  has  to  be  calibrated  with  monoenergetic  beams. 
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At  present  the  experimental  results  were  recorded  by  taking  Polaroid  pictures  from  the 
screen  of  th£  oscilloscopes.  The  picture  quality  is  very  poor  as  seen  from  Figs.  2  through 
4,  which  makes  it  difficult  to  analyze.  The  use  of  film  turned  out  to  be  very  expensive. 
We  urgently  need  to  set  up  the  computer  data  acquisition  system.  The  hardware  interface 
which  was  obtained  without  cost  was  installed  and  software  has  to  be  implemented  and 
tested. 

It  is  not  possible  to  discuss  all  the  diagnostic  system  and  measurements  and  the  data 
acquisition  system.  To  build,  test,  calibrate  and  identify  corrrect  measurements  is  a  major 
part  of  the  scientific  effort. 

The  first  stage  experiment  will  be  carried  out  with  a  mirror  field  of  6  kG,  which  is  about 
half  of  the  design  value  and  just  above  the  ICRF  resonance  magnetic  field.  The  mirror 
coils  can  be  run  at  room  temperature  with  intermittant  cooling.  To  run  at  higher  fields 
the  coils  have  to  be  cooled  down  to  liquid  nitrogen  temperature.  The  rf  power  injection 
will  be  started  at  a  low  power  level  (10-20  kW)  and  short  pulse  length  (10-20  ms).  The 
power  level  and  pulse  length  will  be  increased  gradually  after  the  plasma  properties  are 
fully  understood  and  optimal  operating  conditions  are  obtained.  Study  will  be  made  on 
the  effect  on  heating  with  varying  gas  flow  rates,  field  intensities,  and  microwave  radiation 
feedthrough  locations. 

During  the  course  of  the  heating  experiment  we  will  try  to  study  the  diamagnetic  loop 
for  plasma  3  measurement  and  to  complete  implementation  of  a  computer  data  acquisition 
and  reduction  system.  A  microwave  interferometer  will  also  be  installed  for  plasma  density 
measurement. 

More  diagnostic  equipment  is  needed.  It  typically  takes  at  least  two  months  to  make 
a  given  diagnostic  operational  and  to  understand  the  measurements.  In  addition  to  the 
experimental  study  the  following  specific  tasks  are  planned:  (1)  to  study  the  diamagnetic 
loop  measurements;  the  loop  is  outside  the  chamber  and  may  need  to  be  moved  inside 
the  chamber  (2)  to  acquire  a  spectrometer  to  do  Doppler  broadening  measurements  for 
obtaining  ion  temperature;  (3)  to  construct  and  install  a  second  antenna  in  the  central 
cell  ?nd  to  tune  the  rf  transmitter  for  injecting  the  power  into  the  central  cell;  (4)  to 
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carry  out  heating  experiments  in  both  the  central  cell  and  the  anchor;  (5)  to  make  the 
third  rf  transmitter  operational  for  heating  the  other  anchor;  and  (6)  to  convert  the  pulsed 
microwave  transmitter  to  full  wave. 


7  Instrumentation  Requirements  (1990) 

To  understand  the  plasma  properties  correctly  we  need  many  diagnostic  instruments. 
At  present  we  urgently  need  the  following 

1.  Langmuir  probes 

Langmuir  probes  will  be  used  to  measure  the  plasma  edge  density  and  temper¬ 
ature  at  various  locations.  A  known  profile  will  be  used  to  project  the  center 
conditions.  We  already  have  made  one  set  of  working  probes.  We  need  to  build 
more  probes  and  a  power  supply  to  drive  the  probes. 

2.  Diamagnetic  loop 

A  diamagnetic  loop  will  be  used  to  measure  the  /3  in  the  central  cell  and  one  of 
the  plugs.  3  is  the  ratio  of  plasma  pressure  to  magnetic  field  pressure  P/(B2 / 2p). 
Therefore  we  can  determine  the  over  plasma  pressure  P  =  j3{B2 /2^). 

3.  B  probe  array 

A  B  probe  array  will  be  used  to  monitor  magnetic  fluctuations  and  propagation 
of  rf  waves.  From  this  we  will  study  the  method  of  improving  the  efficiency. 

4.  Bolometer 

A  bolometer  and  electronics  will  be  used  to  monitor  radiation  losses. 

5.  Grid  energy  analyzer 

A  grid  energy  analyzer  and  electronics  will  be  used  to  measure  plasma  potential 
in  the  plug  and  particle  energy  at  the  exhaust  from  which  we  can  determine  the 
exhaust  velocity. 


6.  Visible  light  spectrometer 


A  visible  light  spectrometer  will  be  used  to  measure  ion  temperature  from  the 
Doppler  broadening  of  the  light. 

7.  Data  acquisition  system 

We  obtained  an  agreement  to  share  the  data  acquisition  VAX  computer  with 
the  fusion  energy  research  groups.  The  interface  hardward  and  software  have 
been  installed.  We  are  testing  the  operation.  However,  our  CAM  AC  crater  and 
digitizer  are  too  slow.  High  speed  CAMAC  craters  and  digitizers  are  needed. 

8.  Emissive  probes  to  measure  plasma  potential. 

8  Available  Facility 

The  propulsion  experimental  facility  has  been  set  up  and  is  operational.  The  major 
equipment  is  listed  in  Table  3. 
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Table  3 


Major  Equipment 


ITEM 

RATING 

QUANTITY 

Tandem  Mirror  Device 

Turbo  Pump 

400  liter/sec 

2 

Microwave 

2.4  MHz,  2kW 

1 

Power  Supply 

7,000  kW,  DC,  5000  amp 
(5  MW,  2  sec  pulse) 

2 

18  kW,  DC,  600  amp 

1 

150  kW,  DC,  5000  amp 

1 

3  kW,  DC,  200  amp 

1 

RF  Transmitter 

300  kW,  DC,  2-30  MHz 

3 

Dye  Laser 

1 

Ruby  Laser 

1 

CAMAC  Interface 

1 

Large  Vacuum  Tank 
with  Pumping  System 

1 
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